Field-induced suppression of the 7r-band 
superconductivity and magnetic hysteresis in the 
microwave surface resistance of MgB2 at 
temperatures near Tc 

M. Bonura, A. Agliolo Gallitto and M. Li Vigni 

CNISM and Dipartimento di Scienze Fisiche e Astronomiche, Universita di Palermo, 
Via Archirafi 36, 1-90123 Palermo (Italy) 

G. A. Ummarino 

CNISM and Dipartimento di Fisica, Politecnico di Torino, Corso Duca degli Abruzzi 
24, 1-10129 Torino, Italy 

Abstract. We report on the magnetic- field-induced variations of the microwave 
surface resistance, in a polycrystalline MgB2 sample, at different values of 

temperature. We have detected a magnetic hysteresis in Rs, which exhibits an 
unexpected plateau on decreasing the DC magnetic field below a certain value. In 
particular, at temperatures near Tc the hysteresis manifests itself only through the 
presence of the plateau. Although we do not quantitatively justify the anomalous shape 
of the magnetic hysteresis, we show that the results obtained in the reversible region 
of the Rs{H) curve can be quite well accounted for by supposing that, in this range of 
magnetic field, the tt gap is almost suppressed by the applied field and, consequently, all 
the TT-band charge carriers are quasiparticles. On this hypothesis, we have calculated 
Rs{H) supposing that fluxons assume a conventional (single core) structure and the 
flux dynamics can be described in the framework of conventional models. From the 
fitting of the experimental results, we determine the values of -ff^2(^) temperatures 
near Tc. In our opinion, the most important result of our investigation is that, at least 
at temperatures near Tc, the value of the applied field that separates the reversible 
and irreversible regions of the Rs{H) curves is just i?^2(^); ^ qualitative discussion of 
the possible reason of this finding is given. 



PACS numbers: 74.25.Ha; 74.25.Nf; 74.60.Ge 



Magnetic hysteresis in the niw surface resistance of MgB2 at fields lower than H^2 2 
1. Introduction 

A suitable metliod to investigate fiuxon dynamics in type-II superconductors consists 
in measuring the magnetic-field-induced variations of the microwave (mw) surface 
resistance, Rg [H El El H]. Indeed, the field dependence of Rg in superconductors in 
the mixed state is determined by the presence of fluxons, which bring along normal 
fluid in their cores, as well as the fiuxon motion. The experimental results are generally 
discussed in the framework of the two-fluid model, including the field dependence of the 
quasiparticle density and the effects of the fiuxon motion [5l|6l[7j. 

Studies reported in the literature on the field-induced variations of Rs in MgB2 
have highlighted several anomalies, especially at applied magnetic fields much lower 
than the upper critical field, among which unusually enhanced field-induced mw 
losses [HI El [ini [H] and a magnetic hysteresis of unconventional shape [HI [121 [13]. 
It has been suggested that these anomalies are strictly related to the peculiarities of the 
fiuxon lattice in MgB2. On the other hand, it is by now accepted that fluxons in MgB2 
have a composite structure, being constituted by two concentric cores, one of radius 
(small core), associated with the a gap (Ao-), and the other of radius ^-^ (giant core), 
associated with the tt gap (A^r) [HI [151 [16]. Moreover, the gaps depend very differently 
on the magnetic field: while the a gap closes at the macroscopic upper critical field, Hc2, 
the TT gap is almost closed at magnetic- field values much lower than Hc2- This crossover 
field, at which the contribution to the superconductivity due to the tt superfiuid is 
negligible, has been highhghted in several experiments [HI [TTl [HI [19] and is indicated 
as H^2- Because of the different magnetic- field dependence of the two gaps, on varying 
the field, the structure of the vortex lattice evolves in an unusual way. At low magnetic 
fields, quasiparticles from vr and a bands are trapped within the vortex core, even if on 
different spatial scales. On increasing the field, though cr-band quasiparticles remain 
localized in the small core, giant cores start to overlap because of the field-induced 
suppression of A^; eventually, when is reached, 7r-band quasiparticles are spread 
over the whole sample On further increasing the field, the vr-quasiparticle density 
remains almost unchanged, while the cr-quasiparticle density continues to increase up 
to the macroscopic Hc2 [20j. Only at applied fields higher than H^2y fluxons assume a 
more conventional shape (single core), but they are surrounded by both the condensed 
fluid of the a band and the normal fluid of the vr band. This field-induced evolution of 
the vortex structure is expected to affect both the vortex- vortex and the vortex-pinning 
interactions, making the standard models inadequate to describe the fiuxon dynamics in 
MgB2 in a wide range of magnetic fields. Sarti et al. [21], investigating the mw surface 
impedance of MgB2 films, have shown that at low fields, when the contribution of the 
TT-band superfiuid cannot be neglected, the magnetic-field dependence of the real and 
imaginary components of the surface impedance exhibit several anomalies. However, 
they have shown that for fields higher than a threshold value, lower than Hc2{T), the 
experimental results can be justified in the framework of a generalized two-fluid model in 
which, under the hypothesis that the vr-quasiparticle density has reached the saturation 
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value, the contribution of quasiparticles coming from the vr band is kept constant. 

In this paper, we report on the magnetic-field dependence of the mw surface 
resistance of a polycrystalline sample of MgB2. The field- induced variations of Rs 
have been investigated by the cavity perturbation technique, using a copper cavity, 
at increasing and decreasing the DC magnetic field, Hq. The Rs{Ho) curves exhibit 
several anomalies that cannot be justified in the framework of the standard theories for 
the fiuxon dynamics, among which a magnetic hysteresis of unconventional shape: in the 
decreasing-field branch of the Rs{Hq) curve, we have detected an unexpected plateau 
extending from a certain value of the magnetic field down to zero. The hysteresis is 
detectable up to temperatures close to Tc (T/Tc ~ 0.95); however, the extension of the 
plateau depends on T and, in particular, for T > 0.77 Tc the hysteresis manifests itself 
only through the presence of the plateau. In a previous paper [11], we have extensively 
discussed the anomalies of the Rs{Ho) curves in MgB2 obtained at low temperatures and 
we have ascribed them to the unusual vortex structure in this compound. In this paper, 
we devote attention to the range of fields at which the superconductivity coming from 
the 71 band is almost suppressed. In this region {Hq > expect that the flux 

lines assume a conventional single-core structure and that all the charge carriers coming 
from the vr band are quasiparticles. On this hypothesis, we have modified the expression 
of the complex penetration depth of the mw field, considering that the contribution of 
the TT band to the field-induced energy losses is merely due to the presence of the vr 
quasiparticles and that of the a band is due to both the a quasiparticles and the fiuxon 
motion. We will show that the experimental results can be quite well justified, in the 
framework of the proposed model, for fields greater than a threshold field, which we infer 
to be The fitting of the experimental data allows determining with good accuracy 
the temperature dependence of if^g temperatures near T^, where other techniques fail 
in giving accurate results. Furthermore, our results show that, at least at temperatures 
near Tc, the values of i7c2(^) coincide with the values of the applied magnetic field below 
which the decreasing-field branch of the Rs{Hq) curve exhibits the plateau. 

2. Experimental apparatus and sample 

The field-induced variations of the mw surface resistance have been investigated in a bulk 
sample of Mg^°B2 prepared at the INFM-LAMIA/CNR laboratory in Genova using the 
so-called one-step method [22]. The sample has parallelepiped shape with approximate 
dimensions 2 x 3 x 0.5 mm^; it undergoes a sharp superconducting transition with onset 
Tc ~ 38.9 K and ATc ~ 0.3 K. Figure [1] shows the SEM micrograph of the characteristic 
morphology of the sample surface; it highlights pores of the order of ~ 1 /zm. 

The mw surface resistance has been measured by the cavity-perturbation technique. 
A copper cavity, of cylindrical shape with golden-plated walls, is tuned in the TEqh mode 
resonating at uj/27i ^ 9.6 GHz [Q ^ 40000 at LHe temperature). The sample is located 
in the center of the cavity, by a sapphire rod, where the mw magnetic field is maximum. 
The cavity is placed between the poles of an electromagnet which generates DC magnetic 




fields up to 1 T. Two additional coils, independently fed, allow compensating the residual 
field and working at low magnetic fields. A LHe cryostat and a temperature controller 
allow working either at fixed temperatures or at temperature varying with a constant 
rate. 

The sample and the field geometries are schematically shown in Figure [2^; the DC 
magnetic field, Hq, is perpendicular to the mw magnetic field, H^^. When the sample is 
in the mixed state, the induced mw current causes a tilt motion of the vortex lattice [6]; 
Figure [Sb schematically shows the motion of a flux line. 

(b) 
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t 

Figure 2. (a) Field and current geometry at the sample surface; w « 2 mm,t w 
0.5 mm, h K 3 mm. (b) Schematic representation of the motion of a flux line. 




The surface resistance of the sample is given by 

where Ql is the quality factor of the cavity loaded with the sample, Qu that of the 
empty cavity and F the geometry factor of the sample. 

The quality factor of the cavity has been measured by an hp-8719D Network Analyzer. 
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All the measurements have been performed at very low input power; the estimated 
amplitude of the mw magnetic field in the region in which the sample is located is of 
the order of 0.1 /iT. 

3. Experimental results 

The field-induced variations of Rs have been investigated for different values of the 
temperature. For each measurement, the sample was zero-field cooled (ZFC) down to 
the desired temperature; the DC magnetic field was increased up to a certain value and, 
successively, decreased down to zero. 

Figure [3] shows the field-induced variations of RgsXT = 4.2 K, obtained by sweeping 
Hq from zero to 1 T and back. In the figure, ARs{Ho) = Rs{Ho,T) — R^es and 
/\]^max ^ _ /^^^^^ where Rres is the residual mw surface resistance at T = 2.5 K 
and Hq = 0, and i?„ is the normal-state surface resistance at T = Tc. The inset shows 
the results of the increasing-field branch in a logarithmic scale, which allows identifying 
the value of the applied magnetic field at which Rg deviates from its zero-field value; this 
value should be the first-penetration field of vortices. Hp. The decreasing-field branch 
is characterized by two characteristic fields: i) Hirr, which indicates the value of the 
applied field that separates the reversible and irreversible part of the Rs{Hq) curve; ii) 
H', which identifies the beginning of a plateau in the curve. Magnetic hysteresis in 
Rs is expected as a consequence of the critical state of the vortex lattice; it has been 
detected in other superconductors [H [231 121] and has been ascribed to the different B 
value at increasing and decreasing fields. However, as we will discuss in more detail in 
the next section, the presence of the plateau is puzzling because it should indicate that 
the trapped flux does not change anymore on decreasing the field below H'. 
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Figure 3. Field-induced variations of i?s at T = 4.2 K obtained by sweeping Ho from 
to 1 T (circles) and back (triangles). The inset shows the results in a logarithmic 
scale to better identify Hp. 

For T < 0.77 T^, the decreasing-field branch of the Rs{Hq) curve shows peculiarities 
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similar to those of Figure [31 with the two characteristic fields, though they decrease 
on increasing T. For T > 0.77 T^, Hi„.{T) coincides with H'(T), i.e. the hysteresis 
manifests itself only through the presence of the plateau below H'{T); Figures H] and [5] 
show the field-induced variations of Rs, just in this range of temperatures. In all the 
figures, ARs{Hq) = Rs{Hq,T) — Rres and Ai?^"^ = Rn — Rres- The continuous lines 
reported in the figures are the best-fit curves obtained as described in Section 14.21 

I I I I I I I I I I I I I I I I I I I I 
: 7=30 K 
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Figure 4. Field-induced variations of Rs obtained by sweeping Hq from to 1 T 
(circles) and back (triangles). The continuous lines are the best- fit curves obtained as 
described in Sec. 14.21 

The zero- field value of the increasing-field branch of the Rs{Hq) curve is determined 
by the contribution of quasiparticles due to thermal breaking of Cooper pairs. At 
Hq = Hp, vortices start to penetrate the sample giving rise to an increase of Rg. When 
the applied magnetic field reaches the value of the upper critical field, Rs assumes 
its normal-state value, Figure [6] shows the temperature dependence of Hp and 

Hc2 deduced from the isothermal Rs{Ho) curves. The values indicated in panel (b) as 
full circles have been deduced measuring the magnetic field at which Rg reaches Rn, 
those indicated as open circles have been obtained by fitting the experimental data as 
discussed in Sec. 14.21 In any case, since our sample is a polycrystal the upper critical 
field so determined coincides with H;!f. Dashed line is the curve calculated by the 
generalized two-band Eliashberg theory, as reported in the Appendix. 
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Figure 5. Field- induced variations of Rg obtained by sweeping Hq from to Hmax 
(circles) and back (triangles), with H„iax > Hc2{T). The continuous lines are the 
best-fit curves obtained as described in Sec. 14.21 



From Figure [6^, one can see that Hp exhibits a hnear temperature dependence 
down to low temperatures, consistently with results in MgB2 samples reported by 
different authors [331 EH [35]. As it is known, the first-penetration field may differ 
from Hci if demagnetization and/or surface-barrier effects come into play; however, our 
results are consistent with the lower critical field for MgB2 bulk samples reported in the 
literature [35] . 

By looking at Figures H] and [5l one can see that the hysteretic behaviour of the 
Rs{Hq) curve is detectable up to temperatures very close to Tc, although restricted in a 
narrow field range; only at T > 37.5 K, the Rs{Ho) curve is reversible in the whole range 
of fields investigated. This result suggests that pinning is effective up to temperatures 
close to Tc. 

4. Discussion 

Models for the electrodynamics of superconductors in the mixed state have been 
proposed by several authors [3], HI El El [7], who discuss the field-induced mw losses 
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Figure 6. Temperature dependence of the first-penetration field (a) and the upper 
critical field (b), deduced from the isothermal Rs{Hq) curves. Full circles have been 
directly deduced from the experimental data; open circles have been determined by 
fitting the data. In any case, they identify H^{T). Dashed line is the curve calculated 
by the generalized two-band Eliashberg theory, as reported in the Appendix. 



in different approximations. Coffey and Clem (CC) have elaborated a comprehensive 
theory, in the framework of the two-fluid model of superconductivity, in which, besides 
the effects of the fluxons motion, the field dependence of the densities of the normal and 
condensed fluids are considered [5]. The CC theory has been developed under the basic 
assumption that vortices generate a magnetic induction field, S, uniform in the sample. 
On this hypothesis, the mw surface impedance in the linear approximation {H^^ ^ Hq) 
is given by 

Rs = -/^o^^ Im[A(u;, B, T)] , (1) 

with 



T,, , „ y{B,T) + {i/2)PJuj,B,T) 

^(-.^■r) = yi_2,A^(B.T)AJ,(.,B,T)- 

In Eq. ([2]), 6y{uj, B,T) is the effective complex skin depth arising from the vortex 
motion [51 [25] , while the temperature and field dependencies of the two fluids are taken 
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into account through X{B,T) and 6nf{uj,B,T), given by 



\{B,T) = —=^=, 



(3) 



Snf{uj,B,T) 



2 



(4) 



/ioc^cr„(5,T) ' 



where Aq is the London penetration depth at T = 0; 6o is the normal skin depth at 
T = Tc] ns{B,T) and nn{B,T) are the condensed- and normal-fluid fractions, respec- 
tively; an{B,T) is the normal-fluid conductivity in the superconducting state. 

Equations ([IHl]) have been obtained under the hypothesis that B is uniform inside 
the sample. When the fluxon lattice is in the critical state, the assumption of uniform B 
is no longer valid; as a consequence, the hysteresis in the Rs{Hq) curve cannot be justified 
by Eqs. (1-4). In our field geometry (see Fig. [2^), the effects of the non-uniform B 
distribution on Rg are particularly enhanced because in the two sample surfaces normal 
to the external magnetic field the mw current and fields penetrate along the fiuxon axis 
and, consequently, the mw losses involve the whole vortex lattice. However, in this case, 
one can easily take into account the non-uniform B distribution by calculating a proper 
averaged value of Rg over the whole sample as follows [31 H] 



where S is the sample surface, S is its area and r identifies the surface element. 
It is worth noting that, in order to use Eq. ([5]) for taking into due account the critical- 
state effects, it is essential to know the B profile inside the sample, determined by the 
field dependence of the critical current density. 

Recently, using this method, we have investigated the effects of the critical state 
on the field-induced variation of Rs, at increasing and decreasing fields [31 Hj. We have 
shown that the parameter that mainly determines the peculiarities of the Rs{Hf)) curve 
is the full penetration field, H*. In particular, H* strongly affects the amplitude of 
the hysteresis loop and the shape of the increasing- field branch of the Rs{Ho) curve. 
However, independently of the H* value, the shape of the decreasing-field branch, 
being strictly related to the shape of the magnetization curve, should exhibit a negative 
concavity, with a monotonic reduction of Rs from Hq = Hmax down to Hd- 

The Rs{Hq) curves we have obtained cannot be justified in the framework of 
the model above mentioned, for several reasons. From Figure [31 one can see that 
the application of a magnetic field of ~ IT, which is of the order of ifc2/10, 
causes a Rg variation larger than 30% of the maximum variation. These field-induced 
variations of Rg are much greater than those expected from the models reported in 
the literature [31 IH El El [7| and detected in other superconductors [H [25l [26]. Another 
anomaly concerns the shape of the magnetic hysteresis and, in particular, the presence 
of the plateau in the decreasing-field branch. As above mentioned, the decreasing-field 




(5) 



Magnetic hysteresis in the niw surface resistance of MgB2 at fields lower than H^2 10 

branch of the Rs{Hq) curve is expected to exhibit a monotonic reduction of Rg down to 
Hci- The presence of the plateau is puzzhng; indeed, it would suggest that the trapped 
flux does not change anymore on decreasing the field below H', although this value is 
much larger than the first penetration field. 

We would like to remark that results similar to those reported in this paper have 
been observed in other MgB2 bulk samples, produced by different methods |TT|,[T3]. We 
have extensively discussed the anomalies of the Rs{Ho) curve in Refs. [11], [13], and we 
have suggested that they are ascribable to the unusual properties of the fiuxon lattice 
in two-gap MgB2. This hypothesis has been corroborated by the fact that the results 
obtained in a strongly irradiated MgB2 sample, in which the two gaps merged into a 
single value, have been quite well accounted for in the framework of the model above 
discussed, using Eqs. ([THS]) [HI ETj. In two-gap MgB2 samples, according to Ref. [8], 
the unusually enhanced field-induced variations of Rs at applied fields much lower than 
Hc2 have been qualitatively ascribed to the presence and motion of the giant cores due 
to the TT-band quasiparticles. On the contrary, the origin of the anomalous shape of 
the Rs{Ho) curve is so far not understood. Only in finely powdered MgB2 samples the 
hysteresis assumes a more conventional shape (without plateau), suggesting that the 
presence of the plateau is related to a bulk process [T3] . 

By looking at Figures H] and [5l one can see that the hysteretic behaviour of the 
Rs{Ho) curve is detectable up to T ^ 0.95 Tc, suggesting that pinning is effective up to 
temperatures very near T^. However, on increasing the temperature the irreversibility 
is restricted in a narrow range of fields and, for T > 0.77 Tc, it manifests itself by the 
mere presence of the plateau; in the following, we will focus on the analysis of the results 
obtained just in this temperature range. 

4.I. The model 

The study of the field-induced variations of Rg in MgB2 is made particularly complex 
by the unconventional structure of vortices, which is expected to affect both the vortex- 
vortex and the vortex-defect interactions. However, for Hq > H^2y vortices are expected 
to assume a "single-core" structure; indeed, giant cores overlap and a quasiparticles 
remain localized within the small cores, which will overlap at Hq = Hc2- Starting from 
this picture, we have modified the expression of the complex penetration depth of the 
mw field; we have considered that the contribution of the vr band to the field-induced 
energy losses is merely due to the presence of the n quasiparticles and that of the a 
band is due to both the a quasiparticles and the fiuxon motion. We can neglect the 
effects of the critical state if we look at the reversible region of the Rs{Ho) curve. 

Since in MgB2 the superconductivity is due to the charge carriers coming from the 
two different bands, one can write: Ug = n^ + nl and rin = n^+nf^. For < Hq < Hc2, 
the superelectron fraction, n^, reduces to because = 0. So, \{B,T) is determined 
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by the a band and Eq. can be rewritten as 

X(B T) = = ^0 (Q\ 

y^^WTr) y/n'^AO,T)[l-B/B,,{T)] ' 

In MgB2, due to the presence of normal fluid coming from the two bands, the overall 
normal-state conductivity can be reasonably considered as the sum of the normal-state 
conductivities of the vr and a bands: an = cr^ + cr^. Consequently, the normal skin depth 
can be written as 



So=\ = W 7 r- (7) 



Sarti et al. [2T] have suggested that, since the two bands interact very weakly each other 
[28j, when superconductivity arising from the a band is suppressed the conductivity of 
the quasiparticle fraction can be written as 

T) = < + T) = < + [1 - < (5, T)]< . (8) 

It follows that, in the range of fields considered, Eq. (jl]) can be modified as 

In order to calculate the complex penetration depth of the mw field, besides X{B,T) 
and 6nf{uj, B,T), it is necessary to determine the effective complex penetration depth 
due to the vortex motion (see Eq. [2]). S^{uj,B,T) depends on the relative magnitude 
of the viscous and restoring-pinning forces through the depinning frequency [29] . Since 
we are analyzing the results in a restrict range of temperature near Tc, and considering 
the values of the depinning frequency reported in the literature for MgB2 [101 ED], it 
is reasonable to assume that the vortex motion at our working frequency is ruled by 
the viscous drag force, i.e. vortices move in the fiux-fiow regime. The hypothesis is 
strengthened by the consideration that for Hq > H^2 vortices are surrounded by cr-band 
condensed fluid and vr-band normal fluid; this should reduce the stabilization energy 
and, consequently, the pinning efficacy. In the flux-flow regime, the expression of 5„ 
reduces to 



6,iu;,B,T) = 5oVB/B,2iT) . (10) 

The equations reported up to now do not explicitly consider the anisotropy 
properties of MgB2. In order to take into account the anisotropy, one can assume that 
the polycrystalline sample is constituted by grains with the c-axis randomly oriented 
and suppose that at fixed temperatures the anisotropic Ginzburg-Landau theory can be 
used. By indicating with 9 the angle between Hq and the c-axis of the generic crystallite, 
the upper critical field is given by 

cos2(e) + sin2(e) 
where 7 is the anisotropy factor of the upper critical field. 
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Also the anisotropy of the penetration depth plays a role because it determines the 
sample surface layers in which the mw energy losses occur. Indicating with a the angle 
between H^^ and c one can write 

A^(a) = (A^),,^cos2(a)+72sin2(a), 

where 

It is worth noting that, since in our field geometry Hq ± H^, a ^ 6. Supposing 
Hq = z and H^^ = x, one can easily find a = arccos(sin6'cos(y9), being 9 and cp the polar 
and azimuthal angles. The mw surface resistance of the generic crystallite depends on 
the angles 9 and a between its c-axis and the directions of Hq and if^, respectively; the 
overall surface impedance of the sample can be obtained by integrating over the whole 
solid angle. 

4-2. Analysis of the results 

The expected results of the normalized Rs{Hq) curves depend on {XQ)ab/So, 
ng{0,T), a^/an, 7 and 7a. Most of these quantities are known and/or deducible from 
the experimental data. 

The anisotropy of Ag can be determined using the results of first-principles band- 
structure calculations [36] ; in particular, from the values of the plasma frequencies 
reported by Brinkman et al. [37], one obtains 7a = 6.1. The temperature dependence of 
n" at zero magnetic field has been obtained by using the generalized two-band Eliashberg 
theory, which has already been used with success to study the MgB2 system [38] . 
Figure [7] shows the calculated temperature dependence of n"] details of the calculation 
procedure are reported in the Appendix. In the analysis of our experimental results, we 
have used for n^(0,T) the values reported in Figure [7] and for the anisotropy of Aq the 
value 7a = 6.1. 

The anisotropy of the upper critical field, 7, has been determined both 
experimentally and theoretically [211 [321 13S]; at temperatures near Tc, it should be 
in the range 2 3. In the temperature range 34 -h 37.5 K, H^{T) has been deduced 
from the experimental data (full circles in Fig. [6l); for T < 34 K it has to be taken as a 
parameter. Moreover, we have found that for T > 33 K the value of cr^/cr„ affects very 
little the expected curves; this finding can be ascribed to the fact that at temperatures 
close to Tc, 5nf{uj, B, T) ^ (Jo- 
in order to fit the results obtained in the temperature range 34 -^ 37.5 K, we have 
used for the values reported in Figure [6] as full circles, letting them vary within 
the experimental uncertainty; have let 7 vary in the range 2 3; have taken {XQ)ab/5o 
as free parameter. The best-fit curves are shown as continuous lines in Figure |5l the 
values of the best-fit parameters are reported in Table 1. 

For T < 34 K, our apparatus does not allow measuring the upper critical field, 
which has to be taken as parameter as well; on the other hand, we can use for (AQ)afe/(5o 
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Figure 7. Temperature dependence of at zero magnetic field, calculated by the 
generalized two-band Eliashberg theory, using the procedure reported in the Appendix. 
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Table 1. Values of the parameters that best fit the experimental data of Figs. [Hand O 
the values reported without uncertainty have been imposed. For all the temperatures, 
we have used 7a = 6.f and nj(0, T) of Fig. 7. 

the value we have aheady determined by fitting the data at higher temperatures. At 
T = 30 K, also the value of cr^/cr„ affects in a not-negligible extent the expected curve; 
so, the parameters necessary to fit these data are H;^'^ and 7. The best-fit curves 

obtained at T = 30 K and T = 33 K are shown in Figure H] and the parameters are 
listed in Table 1. 

In the framework of the Eliashberg theory, we have calculated the temperature 
dependence of the upper critical field. We have used the linearized gap equations under 
magnetic field, for a superconductor in the clean limit [39]; details of the calculus are 
reported in the Appendix. In this calculation, the only input parameters are the Fermi 
velocities in the two bands. The calculated values of H^{T) are reported, as dashed line 
in FigureEl they have been obtained with w^^j, = 4.4 x 10^ m/s and fj^^^ = 20 x 10^ m/s. 
^Fab equal to that used in the band-theory calculus of Brinkman et al. [37]; on the 
contrary, fj.^^ is larger than that reported by Brinkman et al., as obtained also by other 
authors [39j . 

From the deduced value of {\Q)ab/5oi one can estimate (Aq)^^ if 5q is known. 
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The value of the normal-state mw surface resistance, we have measured at T = T^, 
is Rn ~ 40 mfi; from the estimated Rn we deduce 60 = 2Rn/ fiQUj ^ 1 fim and, 
consequently, {\Q)ab ~ 220 nm. This value of (Ao)ab is larger than the expected value 
calculated for a single crystal; we think that this is due to the porousness of the sample 
and/or the roughness of its surface (see Fig. [1]). Indeed, the surface roughness enlarges 
the effective area of the sample in which the mw energy losses occur, giving rise to an 
overestimated i?„ value, with consequent enlarged value of the deduced 60 and (Aq)^^- 

Although the values of the best-fit parameters have been obtained by fitting the 
results for Hq > H', the expected curves of Figures H] and [5] are reported in the whole 
range of fields investigated. In the framework of our model, the expected curves should 
properly describe the experimental results only for Hq > H^2- -^^ lower fields, it is 
expected that the theoretical curve overestimates Rs{Ho) because our model assumes 
that all the vr-band carriers are quasiparticles. Therefore, the Ho value at which the 
theoretical curve begins to fit the experimental data can be taken as ifjgl^)- Using 
this criterion, we have deduced the temperature dependence of the range of 

temperatures investigated. The results are shown in Figure |8] as open circles; in the same 
figure, we report the values of ifjgl^) obtained by Daghero et al. (full triangles) [H] 
and Samuely et al. (full squares) [19] from point-contact spectroscopy. 




T/r 



Figure 8. Temperature dependence of Open circles are the results we have 

obtained as described in the text; full triangles and full squares are the results obtained 
by Daghero et al. [18] and Samuely et al. 19J, respectively. 

The values we have obtained for H^2(T) match with those reported in Refs. [18] and 
[T9] for lower temperatures; they give information on the temperature dependence of 
H^2 ^ range of temperatures close to Tc, where conventional techniques fail in giving 
accurate results. 

As one can see from Figures H] and O we have obtained good fits for Hq > H'(T); i.e. 
H'^2{T) coincide with the values of Hq at which the decreasing- field branch of the Rs{Hq) 
curves deviates from the increasing-field one, within the experimental uncertainty. This 
finding suggests that the features of the Rs{Hq) curves are strongly related to the 
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suppression of the vr gap, due to the apphed magnetic field. This is a curious result 
that needs further investigation; however, some considerations can be done. 

As already mentioned, results very similar to those reported in this paper have 
been obtained in all the MgB2 bulk samples we have investigated [HI [121 [13], except 
in a strongly irradiated MgB2 bulk sample, in which the two gaps merged into a single 
value [ll]. Moreover, the anomalous plateau in the decreasing-field branch of the Rs{Ho) 
curve has not been observed in samples of MgB2 consisting of fine powder [13]. In order 
to propose a possible explanation for the different results in bulk and powder samples, 
it is worth to recall that, in our field configuration (HqJ^Hi^), the mw current induces a 
tilt motion of the whole vortex lattice (see Fig. [2]); so, our measurements allow detecting 
the response of the fiuxon lattice also in the interior of the sample (far from the edges 
from which the DC magnetic field penetrates). On the other hand, in the framework 
of the critical state, the main property distinguishing the magnetic response of samples 
having small or large dimensions is the residual magnetic induction far from the edges, 
after the applied magnetic field has reached values larger than the full penetration field. 
We infer that the anomalous shape of the hysteresis in Rs{Hq) detected in the MgB2 
bulk samples is related to a process occurring in the interior of the sample, probably 
due to the higher value of the local magnetic field. 

The finding that the value of the applied field that separates the reversible and 
irreversible regions of the Rs{Ho) curves coincides with H^2(T) can be justified by the 
consideration that for Hq > vortices are surrounded by a-band condensed fiuid and 
vr-band normal fiuid; this should reduce the stabilization energy and, consequently, the 
pinning efficacy. The hysteretic behaviour at Hq < H^2 could be justified by supposing 
that, after the sample has been exposed to fields higher than H^2y on reducing Hq, the 
induction field inside the sample, far from the surface, is large enough to destroy the 
superconductivity in the vr band but, contemporaneously, the pinning is effective in the 
sample regions near the surface, where the local magnetic induction is smaller than /zq 
H^2- When the DC field is reduced below H^2y the vortices near the surface experience 
the restoring pinning force which, hindering the free exit of fiuxons, may keep the inner 
induction field at values ~ /xq H^2'^ this may imply high dissipation due to both the high 
quasiparticle density and the free motion of fiuxons. These considerations, though may 
justify the presence of the magnetic hysteresis, do not quantitatively account for the 
plateau; the latter can be justified only if one assumes that the region of the sample 
in which the fiuxon lattice is in the critical state, after it has been exposed to applied 
fields higher enough than H^2^ is much smaller than that in which the local induction 
maintains the value /io -^^2- 

Conclusion 

We have investigated the field-induced variations of the microwave surface resistance 
of a polycrystalline sample of MgB2, at fixed temperatures. The mw surface resistance 
has been measured by the cavity-perturbation technique, at increasing and decreasing 
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the DC magnetic field in the range 1 T. The Rs{Hq) curves exhibit a magnetic 
hysteresis of unconventional shape: in the decreasing- field branch of the Rs{Hq) curve, 
we have detected an unexpected plateau extending from a certain magnetic-field value 
down to zero. The hysteresis is detectable up to temperatures close to Tc {T /T^ ^ 0.95); 
however, for T > 0.77 Tc the hysteresis manifests itself only through the presence of the 
plateau. 

It has been shown by different authors that the standard models for the investigation 
of fluxon dynamics are inadequate to justify the results in MgB2 in wide ranges of 
fields and temperatures, due to the unusual vortex structure (double core). Our results 
confirm this idea; however, in this paper, we have devoted attention to the range of 
fields at which the superconductivity coming from the tt band is almost suppressed. 
In this region {Hq > H^2)^ expect that the flux lines assume a conventional single- 
core structure and all the charge carriers coming from the tt band are quasiparticles. 
On this hypothesis, we have modified the expression of the complex penetration depth 
of the mw field, considering that the contribution of the tt band to the field-induced 
energy losses is merely due to the presence of the tt quasiparticles and that of the a 
band is due to both the cr quasiparticles and the motion of conventional (single core) 
fiuxons. By using the generahzed two-band Eliashberg theory, we have calculated the 
temperature dependence of the a-band quasiparticle density, at zero field, which has 
been successfully used for determining the expected Rs{Ho, T) curves. In the framework 
of the same theory, we have also calculated the temperature dependence of the upper 
critical field and compared it with that deduced from the experimental results, in the 
range of temperatures investigated. 

We have quantitatively analyzed the experimental results obtained at temperatures 
near Tc {T/Tc > 0.77), where our experimental apparatus allows reaching DC magnetic 
fields larger enough than H^2- Since the used model is valid only for Hq > H^2: by 
looking at the field range in which the expected results fit the experimental data, we 
have determined with good accuracy the temperature dependence of H^2 temperatures 
near Tc, where the standard techniques fail in giving accurate results. Our results show 
that the applied magnetic field that separates the reversible and irreversible region of 
the Rs{Ho) curves is just H^2{T)- This interesting result seems to link the magnetic 
field that suppresses the superconductivity of the tt band to the pinning efficacy. We 
have suggested that for Hq > H^2 the Rs{Ho) is reversible because the vortex cores are 
surrounded by both cr-band condensed fluid and 7r-band normal fluid; this reduces the 
stabilization energy and, consequently, the pinning efficacy. The hysteresis observed for 
Hq < Hc2 can be justified assuming that vortex pinning is effective up to temperatures 
near Tc when charge carriers from both vr and a bands contribute to the condensed fluid. 

The reason for which the hysteretic behaviour of the Rs{Hq) curve at temperatures 
near Tc manifests itself only through the presence of the plateau, extending from H^2 
down to zero, is not fully understood and further experimental as well as theoretical 
investigation is necessary. A possible explanation is that, after the sample has been 
exposed to DC magnetic flelds higher that 11^2, on reducing Hq down to zero, in a wide 
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inner region of the sample the local induction field maintains values ~ /io ^^"^ only 
in a narrow region near the sample surface, where the local magnetic field is lower than 
/io the critical state develops. 

Appendix A. Calculus of the superfiuid density at zero magnetic field 

The generalization of the Eliashberg theory for systems with two bands has already 
been used with success to study the MgB2 system [38]. To obtain the gaps, and then 
the superfiuid densities n"(T), n1{T) and the critical temperature, within the s-wave 
two-band Eliashberg model, one has to solve four coupled integral equations for the 
gaps Ai{iujn) and the renormalization functions Zi{iu}n), where i = a, tt is a band index 
and LJn are the Matsubara frequencies. Including in the equations the non-magnetic- 
and paramagnetic- impurity scattering rates in the Born approximation, r!j! and r*i, one 
obtains 

UJnZi{iUn) = U;„ + TtT ^ Aij{iLJn - iUm)N^zi^Um) + 

+ Y.^n + r^i)ivi(^^n) (A.i) 

j 

Zi{iuJn)Ai{iuJn) = nT'^[Aij{iuJn - iuj,m) - fJ'iji^^c)] ■ 

m,j 

■ eilco^l - cojNiitcOm) + - rDNiitcOr.) (A.2) 

j 

where 9 is the Heaviside function, Uc is a cut-off energy, and 

r+oo 

Aij{iujn - iujm) = / dujaijF{uj)/[{uJn - iOmf + Uj"^] , 
Jo 

The solution of the Eliashberg equations requires as input: 

i) the four electron-phonon spectral functions afj{uj)F{u!); 

ii) the four elements of the Coulomb pseudopotential matrix ^*{u!c); 

iii) the four non-magnetic- impurity scattering rates F^; 

iv) the four paramagnetic-impurity scattering rates F^. 

In our calculus is always F^j! = Fj' = 0. The four spectral functions afj{u!)F{u), which 
were calculated for pure MgB2 in Ref . [IQ] , have the following electron-phonon coupling 
constant: A^„=1.017, A^^=0.448, A^^=0.213 and A^^=0.156 

As far as the Coulomb pseudopotential is concerned, we use the expression 
calculated for pure MgB2 [H] 

fJTT 
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c2 



2.23 



Nf,{Ep) N%{Ep) 



2.48 



(A.3) 



where fi{uJc) is a free parameter, NlflEp) is the normal density of states at the Fermi 
level of the i-band and Nj^^Ep) = N^{Ef) + A^^(^f)- 

For obtaining the experimental critical temperature we fix fi{uJc) = 0.03153 with cut-off 
energy Uc = 450 meV and maximum energy 500 meV. In all our calculations we use 
N^{Ef) = 0.30 states/(cell eV) and N^{Ef) = 0.41 states/(cell eV) 



Appendix B. Calculus of the temperature dependence of the upper critical 
field 

In order to calculate the upper critical field we have used the linearized gap equations 
under magnetic field, for a superconductor in the clean limit (negligible impurity 
scattering) [39] 

uj„Zi{iuJn) = ujn + nT'^Aij{iuJn - ium)sign{um) (B.l) 

Zi{iuJn)Ai{iuJn) = 7rT^[Aij{iuJn - iujm) - Kj(t^c)] ■ 

■ 0{\uJc\ - U}m)Xj{i^m)Zj{iuJm)^j{iuJm) (B.2) 

XjiiuJm) = —T= / rfgexp(-g^) • 
■ tan 



{iujm) \ + iHBHc2sign{uJrn) 

(B.3) 

with Pj = TT{vp)'^Hc2/ {2^0)1 is the Fermi velocity of band j and ifc2 is the upper 
critical field. 

In these equations, the bare Fermi velocities are the input parameters for calculating 
the upper critical field as a function of temperature. To obtain the theoretical curve 
reported in Figure El we have used Vp^f^ = 4.4 x 10^ m/s and Vp^^^ = 20 x 10^ m/s. 
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